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ABSTRACT 

We present and analyse a catalogue of 233 pulsars with proper motion measurements. 
The sample contains a wide variety of pulsars including recycled objects and those as- 
sociated with globular clusters or supernova remnants. After taking the most precise 
proper motions for those pulsars for which multiple measurements are available, the 
majority of the proper motions (58%) are derived from pulsar timing methods, 41% 
using interferometers and the remaining 1% using optical telescopes. Many of the 1-D 
and 2-D speeds (referring to speeds measured in one coordinate only and the magni- 
tudes of the transverse velocities respectively) derived from these measurements are 
somewhat lower than earlier estimates because of the use of the most recent electron 
density model in determining pulsar distances. The mean 1-D speeds for the normal 
and recycled pulsars are 152(10) and 54(6) kms respectively. The corresponding 
mean 2-D speeds are 246(22) and 87(13) kms" 1 . PSRs B2011+38 and B2224+64 have 
the highest inferred 2-D speeds of ~ 1600 kms" 1 . We study the mean speeds for dif- 
ferent subsamples and find that, in general, they agree with previous results. Applying 
a novel deconvolution technique to the sample of 73 pulsars with characteristic ages 
less than 3Myr, we find the mean 3-D pulsar birth velocity to be 400(40) kms -1 . 
The distribution of velocities is well described by a Maxwellian distribution with 1-D 
rms a — 265 km s . There is no evidence for a bimodal velocity distribution. The 
proper motions for PSRs B1830— 08 and B2334+61 are consistent with their proposed 
associations with the supernova remnants W41 and G114.3+0.3 respectively. 
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1 INTRODUCTION 

Neutron stars are high-velocity objects. From observations 
of radio pulsars, it has long been known (Gunn & Ostriker 
1970) that as a population they are generally moving much 
faster than their presumed progenitor population, the mas- 
sive OB stars. While the astrophysical applications of a large 
sample of pulsar velocities are manyfold, most importantly 
the physical mechanism for the high velocities is not well un- 
derstood. Lai, Chernoff & Cordes (2001) have summarised 
the various discussed causes for high pulsar velocities. These 
mechanisms, in general, only predict a limited range of ve- 
locities, so that a large sample of pulsar velocities can in 
principle be used to test the mechanisms. 

Pulsar velocities are determined from measurements of 
their proper motion and distance. Prior to the recent pub- 
lication by Hobbs et al. (2004; hereafter Paper I) approxi- 
mately 138 pulsar proper motions had been measured, with 
68 having values greater than 2a. Most millisecond pulsar 
proper motions had been obtained using timing methods 
(e.g. Nice & Taylor 1995; Toscano et al. 1999b; Wolszczan 
et al. 2000 and Freire et al. 2003) . The poor rotational stabil- 



ity observed for young pulsars 1 (referred to as timing noise) 
affected early measurements of proper motions using tim- 
ing methods; for these, interferometers had been used to 
obtain the proper motions (for example, Lyne, Anderson & 
Salter 1982; Bailes et al. 1990; Harrison, Lyne & Anderson 
1993; Fomalont et al. 1997 and Brisken et al. 2002). In Pa- 
per I we described a new method, based on the modelling 
and removal of timing noise from the timing residuals of 
young pulsars using harmonically related sinusoids, that al- 
lowed accurate proper motions to be determined using stan- 
dard timing methods. These proper motions were shown to 
be consistent with proper motion measurements obtained 
using interferometers. We obtained proper motions for 301 
pulsars with 164 having values greater than 2a or with an 



We use the words 'young' and 'old' throughout this paper to 
refer to pulsars with small and large characteristic ages, r c = 
P/2P. The characteristic age assumes that the current rotational 
period of the pulsar is much greater than the initial period and 
that pulsar slow-down is due to a constant magnetic dipole field. 
The characteristic age may not be a good estimation of the true 
age if dipole braking is not strictly followed. 



2 G. Hobbs et al. 




Figure 1. The Galactic motions of the pulsars in our sample. A pulsar is currently at the position indicated by a circle and the track 
indicates its motion for the last 1 Myr assuming no radial velocity. 



uncertainty less than 15masyr _1 in at least one coordinate. 
For 87 of these pulsars, these provide more precise proper 
motions than the earlier measurements. 

Direct distance measurements have only been obtained 
for a small number of pulsars. The majority of distances 
are estimated from the dispersion measure and a Galactic 
electron density model. Using the Taylor & Cordes (1993; 
hereafter TC93) model, Lyne & Lorimer (1994) found the 
mean pulsar birth velocity 2 to be 450(90) kms -1 . Recently, 
Cordes & Lazio (2002; hereafter CL02) provided an updated 
model which, on average, predicts somewhat smaller dis- 
tances than TC93 which will clearly have an impact on the 
calculated velocities. Hereafter, we designate the velocities 
derived from the two models as V TC and V CL . 

Clearly, with these new proper motion determinations 
and a new electron density model it is productive to revisit 
the statistics of pulsar velocities. In Section 2, we describe 
the sample of proper motions used in this work which com- 
bines new results published in Paper I with other proper 
motion values in the literature. In Section 3, we highlight 
the effect of using the CL02 electron density model to ob- 
tain pulsar distances and velocities. In Section 4 we examine 
the statistical properties of various sub-samples of the ob- 
served sample. Following a brief discussion of the motion 
of pulsars in the Galactic plane in Section 5, we estimate 
the 3-D birth speed distribution of non-millisecond pulsars 



2 Throughout this paper, we will present the uncertainties on 
parameters as a value in parenthesis after each quantity. The value 
represents the error (at the 68% confidence level) in the least 
significant digit. 



in Section 6. Finally, in Section 7, we summarise our main 
results and conclusions. 



2 THE PROPER MOTION SAMPLE 

Proper motions were selected from the literature with more 
precise measurements taking precedence if multiple mea- 
surements exist for a specific pulsar. In practice, these proper 
motions were obtained from the ATNF pulsar catalogue 3 
(Manchester et al. 2005) or from Paper I. These proper mo- 
tion measurements all have values greater than 2a or have 
an uncertainty less than 15masyr _1 . As discussed in Pa- 
per I, for a pulsar lying in the ecliptic plane it is not pos- 
sible, using timing methods, to obtain a precise estimate of 
its proper motion in ecliptic latitude. Therefore, only the 
longitudinal component of the proper motion has been mea- 
sured for many of the pulsars in the sample. The timing so- 
lutions for the 87 pulsars obtained from Paper I have been 
updated using the most recent data available from the Jo- 
drell Bank Observatory data archive. The timing solutions 
were obtained in an identical manner to that described in 
Paper I. The most recent observations used in these timing 
solutions are from February 2004. 

The resulting sample of 233 proper motions is provided 
in Table 1. As detailed in the caption, Table 1 also includes 
the 2-D speed (often referred to as the transverse speed), 
Vf L = HtotD CL , and the difference in 2-D speed between 
CL02 and TC93, V^ c - V$ L . 

3 http://www.atnf.csiro.au/research/pulsar/psrcat 
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Table 1. Pulsar proper motions. For each pulsar, the two components of its proper motion are provided followed by a flag stating 
whether these are given in celestial (C) or ecliptic (E) coordinates. The remaining columns provide the total proper motion, a reference 
for the measurement (see the end of the table for the full bibliographic reference), whether the proper motion was determined using 
interferometric (I), timing (T) or optical (O) means, the distance according to the CL02 distance model or from independent distance 
estimates if available, the magnitude of the transverse velocity assuming the CL02 or independent distance estimates (ind.) and the 
difference AVt = ~V^ C — V^ L between the velocities as measured assuming the CL02 and TC93 distances. 



PSR 


Mlong 


Wat 


Coord. 


Mtot 


Ref. 


Obs. 


Dist. 








(mas yr _1 ) 


(mas yr" 1 ) 




(mas yr" 1 ) 




Type 


(kpc) 


(kms- 1 ) 


(kms" 1 ) 


B0011+47 


19.3(18) 


-19.7(15) 


C 


27.6(17) 


1 


I 


1.56 


204 


34 


T3nnoi *70f~ 1 


5.Z(5 ) 


-3.2(4) 




6.1(5) 


Z 


1 


k nn 
5.UU 


1 A c; 
145 


ind. 


15UUZ1 — i ZD 


3.8(4) 


O At A\ 

-2.4(4) 




A c; (A \ 
4.5(4 ) 


Z 


T 
1 


k nn 
5.UU 


1 fl7 

1U / 


ind. 


T3nnoi 70T? 


b.Z(5 ) 


-2.7(5) 




O.o(5 ) 


Z 


1 


k nn 
5.UU 


1 fii 
ibi 


ind. 


Dnnn tot? 
15UUZ1 — 1 Zr 




— 2.9(6) 


o 


5.4(b) 


2 


1 


5.00 


128 


ind. 


B0021-72G 


4.2(14) 


-3.7(12) 


c 


5.6(14) 


2 


T 


5.00 


133 


ind. 


T3nnoi 70U 
15UUZ1 — i Zti 


c; 1 /1 ^ )^ 
5.1(1Z ) 


q fi /1 o"\ 
— o.b(lZ ) 




fi 0/1 Q~\ 

o.z(lz ) 


z 


1 


k nn 
5.UU 


1/17 
14 / 


ind. 


nnnn tot 

BUUzl — 7zl 


4.9(16) 


— 3.9(19) 


C 


6.3(18) 


2 


T 


5.00 


149 


ind. 


T3n001 70T 

15UUZ1 — ( ZJ 


5.o5(lo ) 


Q /I O/l e;^ 
— a.4Z(15 ) 




Q K/l /I \ 

O.o5(14 ) 


Z 


1 


k nn 
5.UU 


1 c;n 
15U 


ind. 


15UUZ1 — / ZiN 


6.8(13) 


-1.1(14) 


o 


6.9(14) 


2 


1 


5.00 


164 


ma. 


J0024-7204O 


4.6(11) 


-2.9(9) 


c 


5.4(11) 


2 


T 


5.00 


128 


ind. 


TfinO/l 700/1 T T 
JUUZ4 — / ZU4U 




-4.5(/) 




fi s/7^ 
o-0(7) 


9 
Z 


1 


^ nn 
5.UU 


1 fii 
ibi 


ind. 


JUUo4— U5o4 


0.5(15 ) 


q 1 / n^ 
-31(9) 


L" 1 

rj 


31(9) 


Q 
O 


1 


n K/i 
U.54 


7n 

/y 


fi k 
b5 


rSUUol — U / 


Q /1 r ^ 
-8(lZ) 




i_r 

Tj 




Q 

o 


1 


U.41 






rSUUo / +OD 


q / i n\ 
3(10) 


7/1 Q\ 

-7(13) 




8(13) 


Q 

o 


1 


q nc; 


lib 


54 


B0114+58 


7(5) 


17(6) 


E 


18(6) 


3 


T 


2.23 


190 


-8 


Tfl"l Q/l OOQ7 


8(d) 


-21(5) 




ZZ\p ) 


Q 
O 


1 


U.5b 


5o 


1 07 
1Z ( 




1 1 

-11(5) 


-19(5) 


r 1 


ZZ(0 ) 


1 


T 

1 


Z.OO 


Qnn 
oUU 


1 


pni no j_ en 


7/1 o^ 

7(12) 




TT 

JZj 




Q 
O 


T 
1 


O 1 o 

z. ly 






T3fl1 /I A i c;o 


— o.o(iy ) 


1 n/ q~\ 
10(3) 




1 1 /Q\ 
11(d) 


Q 

o 


1 


Z.ZZ 


lib 


1 fi 

— lb 


B0149-16 


3.1(12) 


-27(2) 


c 


27(2) 


1 


I 


0.51 


66 


36 


RnonQ /in 




i 5(o5 J 






5 


1 


n c;s 






JUZlo+4ZoZ 




-4(7) 




5(b) 


Q 

o 


1 


Z.b / 


DO 


7ir 

/ 5 


15UZ04 — Oo 




7o/i ^ 

/U(15j 






5 


1 


n 7*5 






nnom i 1 n 


6(7) 


-37(4) 




37(5) 


6 


I 


0.62 


109 


58 


B0320+39 


16(6) 


-30(5) 


C 


34(6) 


4 


I 


1.01 


163 


77 


r5Uozy+54 


17.0(3) 


n c //i > 
— y.5(4) 




19.5(4) 


/ 


T 
1 


1.06 


98 


™ A 

ma. 


n? no q 1 i 
r5Uool+4o 


A /Q\ 

-4(3) 


5(8) 


Hi 


6(7) 


3 


1 


1.64 


47 


13 


"DflQKQ 1 KO 




- y (D) 


T? 

H; 


1 1 f«\ 
11(b) 


Q 
O 


1 


O 7Q 
Z. IO 


1 A K 

145 


i no 
1UZ 


r5Uo55+54 


n on/ 1 q\ 

y.zu(io ) 


8.2(4) 




1 o q/q^ 
lz.o(o ) 


O 

o 


T 
1 


i 1 n 
1.1U 


b4 


( ^ / 


B0402+61 


23(7) 


61(12) 


E 


65(12) 


3 


T 


2.12 


653 


287 


T3ri/i i n i fin 


-4(5) 


1 0/7\ 

-12(7) 




1 Q / 7\ 

13(7) 


Q 
O 


1 


1 1 


/ O 


OA 
z4 


JU4o / — 4 / 10 


i oi /i t;7/n^ 
1Z1.45 / (y J 


71 /1/ln/n^ 

— / i.44y(y ) 




i4u.y i4^y j 


y 


1 


U.14 


y4 


ind. 


B0447— 12 


10(11) 




E 




3 


T 


1.89 






BU40U— lo 


1 o/q\ 
12(8) 








i n 
1U 


T 
1 


Z.OD 






B0450+55 


52(6) 


-17(2) 


c 


55(6) 


4 


I 


0.67 


175 


31 


B0458+46 


"8(3) 


8(5) 


c 


11(5) 


4 


I 


1.39 


72 


20 


B0523+11 


30(7) 


-4(5) 


c 


30(7) 


4 


I 


3.11 


442 


650 


B0525+21 




7(9) 


c 




4 


I 


1.61 






B0531+21 


-13(2) 


7(3) 


c 


15(3) 


11 


O 


2.00 


142 


ind. 


J0538+2817 


-30(6) 




E 




3 


T 


1.23 






B0540+23 


19(7) 


12(8) 


c 


22(8) 


4 


I 


2.06 


215 


154 


B0559-05 


18(8) 


-16(7) 


c 


24(8) 


4 


I 


3.93 


447 


411 


B0609+37 


5(3) 


5(12) 


E 


7(9) 


3 


T 


0.87 


29 


21 


J0613-0200 


2.0(4) 


-7.0(10) 


c 


7.3(10) 


12 


T 


1.71 


59 


17 


B0611+22 


-4(5) 


-3(7) 


c 


5(6) 


4 


I 


2.08 


49 


63 


J0621+1002 


3.5(3) 


-0.3(9) 


c 


3.5(4) 


13 


T 


1.36 


23 


9 


B0621-04 


0(10) 




E 




3 


T 


2.78 






B0626+24 


-7(12) 


2(12) 


c 


7(12) 


4 


I 


2.24 


74 


81 


B0628-28 


-44.6(9) 


20(3) 


c 


48.7(13) 


1 


I 


1.45 


335 


162 



G. Hobbs et al. 



PSR 


Mlong 


Mlat 


Coord. 


Mtot 


Rcf. 


Obs. 


Dist. 




AV T 




(mas yr -1 ) 


(mas yr -1 ) 




(mas yr _1 ) 




Type 


(kpc) 


(kms ) 


(kms -1 ) 


J0633+1746 


138(4) 


97(4) 


C 


169(4) 


14 


O 


0.16 


128 


ind. 


B0643+80 


19(3) 


-1(3) 


C 


19(3) 


4 


J 


1.55 


140 


164 


B0656+14 


44 1 (61 


—2 4(31 


c 


44.2(6) 


15 




0.29 


59 




B0655+64 


-6(5) 


-8(8) 


E 


10(81 


3 


T 


0.49 


23 





1(171 1 —fis^n 


— 15.7(5) 


15.3(6) 




21.9(6) 


12 


T 


0.86 


90 


19 


B0736-40 


-14.0(12) 


-13(3) 


c 


19.0(17) 


1 


I 


2.64 


238 


756 


B0740-28 


-29(2) 


4(2) 


c 


29.3(2) 


10 




2.07 


287 


-25 


10751+18(17 


—0.2(10) 




E 




3 


T 


1.15 






B0751+32 


-4(51 


7(31 


c 


8(41 


4 




1.52 


58 


91 


B0756— 15 


1(4) 


4(6) 


c 


4(6) 


1 




2.95 


56 


15 


B0809+74 


24.02(9) 


-44.0(4) 


c 


50.1(4) 


7 


I 


0.43 


102 


ind. 


B0818-13 


9(9) 


-47(6) 


c 


48(7) 


10 


[ 


1.99 


453 


105 


R0820+02 


5(111 


-1 (8) 

1(0) 




5(111 


4 




1.01 


24 


10 


B0823+26 


61 (31 


—90(2) 




109(3) 


5 




0.36 


186 




B0833— 45 


—49.60(6) 


29.80(10) 


c 


57.86(8) 


16 




0.29 


80 


ind. 


B0834+06 


2(5) 


51(3) 


c 


51(4) 


6 


I 


0.65 


157 


17 


B0835-41 


-2.3(18) 


-18(4) 


c 


18(4) 


1 


[ 


1.60 


137 


225 


Rnqnfi— 1 7 


27(1 1 1 


— 40(1 1 "I 




48(1 1 1 


4 




0.91 


207 


—64 


B0919+06 


18.35(6) 


86.56(12) 


c 


88.48(12) 


17 




1.20 


503 


ind. 


B0940+16 




9(111 


c 




4 




0.81 






B0942-13 


— 


-22(14) 


c 


— 


4 


I 


0.67 


— 


— 


B0943+10 




-21(12) 


c 




6 


j 


0.64 






R0950+08 


— 2 09(81 


29 46 (71 


c 


29 53(81 


7 




0.26 


36 


inZ 


J1012+5307 


2 4(21 


— 25.2(2) 


c 


25.3(2) 


18 


T 


0.41 


49 


13 


11 1122+1 (1(11 


— 16.3(6) 




E 




3 


T 


0.45 






J1024-0719 


-41(2) 


-70(3) 


c 


81(3) 


12 


T 


0.39 


150 


-15 


B1039-19 


-1(3) 


14(5) 


c 


14(5) 


1 


I 


1.48 


98 


113 


J1045— 4509 


-5(21 


6.0(10) 


c 


7.8(15) 


12 


T 


1.96 


72 


47 


B1112+50 


22(3) 


— 51(31 


c 


56(3) 


4 


I 


0.45 


119 


24 


B1114— 41 


-1(5) 




c 




1 


I 


1.47 






B1133+16 


-74.0(4) 


368.1(3) 


c 


375.5(4) 


7 


I 


0.36 


641 


ind. 


B1237+25 


-106.82(17) 


49.92(18) 


c 


117.91(18) 


7 


I 


0.86 


481 


ind. 


B1254— 10 


-6(81 




E 




3 


T 


1.55 






B1257+12 


46 40(101 


—82.2(2) 


c 


94 39(191 


19 


T 


0.77 


345 


ind 


B1309— 12 


-1(61 




E 




3 


T 


3.14 






B1322+83 


-53(20) 


13(7) 


c 


55(20) 


4 


I 


0.76 


198 


3 


B1325-43 


3(7) 


54(23) 


c 


54(23) 


1 


I 


1.34 


343 


243 


B1426— 66 


—31(51 


— 21(31 


c 


37(51 


20 


I 


1.00 


175 


140 


R1 AdQ — fid 




— 21.3(8) 


Q 


26.6(9) 


20 


j 


2.08 


262 


—30 


J1455— 3330 


5(61 


24(12) 


c 


25(12) 


12 


T 


0.53 


63 


25 


B1451-68 


-39.5(4) 


-12.3(3) 


c 


41.4(4) 


20 


I 


0.45 


88 


ind. 


B1508+55 


-70.6(16) 


-68.8(12) 


c 


98.6(15) 


1 


I 


0.99 


463 


444 


J1518+4904 


4.6(6) 


-8.1(6) 


E 


9.3(6) 


3 


T 


0.63 


28 


3 


B1534+12 


1.340(10) 


-25.05(2) 


c 


25.09(2) 


21 


T 


0.95 


113 


-6 


B1541+09 


-7.3(10) 


-4.0(10) 


c 


8.3(10) 


1 


I 


3.49 


137 


-41 


B1540-06 


-17(3) 


-4(3) 


c 


18(3) 


1 


I 


0.72 


61 


37 


B1552-31 


-9(8) 




E 




3 


T 


2.33 






B1556-44 


1(6) 


14(11) 


c 


14(11) 


10 


I 


2.29 


152 


-47 


J1603-2531 


-10(4) 




E 




3 


T 


1.87 






B1600-27 


-21(11) 




E 




3 


T 


1.56 







A statistical study of 233 pulsar proper motions 



PSR 


Wong 


Wat 


Coord. 


Wot 


Rcf. 


Obs. 


Dist. 


yCL 


A V T 




(mas yr -1 ) 


(mas yr ) 




(mas yr ) 




Type 


(kpc) 


(kms -1 ) 


/i 1 \ 

(kms 1 ) 


J1603-7202 


-3.5(3) 


-7.8(5) 


C 


8.5(5) 


12 


T 


1.17 


47 


19 


B1600-49 


-30(7) 


-1(3) 


C 


30(7) 


20 


I 


5.10 


725 


-215 


R1604— on 


-1C141 


-7(91 


c 


7(10) 


g 


1 


0.67 


22 


— 3 


Ri Ron— riQ 


— 1 7f1 01 




E 




3 


T 


>50.00 






B1620— 26 


— 13.4(10) 


—25(5) 


c 


28(5) 


22 


T 


1.80 


239 


ind 


J1640+2224 


-0.1(3) 


-13.3(7) 


c 


13.3(7) 


19 


T 


1.16 


73 


2 


J1643-1224 


3.0(10) 


-8(5) 


c 


9(5) 


12 


T 


2.41 


103 


105 


RI 642—03 


—3.7(15) 


30 OH 61 




30 2(1 61 


I 


I 


1.12 


160 


256 


RI 70(1— 1 8 






E 




3 


T 


1.48 






B1700— 32 


-31(14) 




E 




3 


T 


2.33 






B1702-19 


-66(5) 


— 


E 


— 


3 


T 


0.88 


— 


— 


J1709+2313 


-3.2(7) 


-9.7(9) 


C 


10.2(9) 


23 


T 


1.41 


68 


20 


R1706— 16 

iJl f U "J _LU 


3(9) 


0(14) 


c 


3(91 


10 


I 


0.83 


12 


6 


11 71 3+0747 


4 Ol^l 


— 4 1 (1 01 


Q 


6 4(7) 


24 


T 


1.11 


34 






— 1 C41 


— 26(5) 


Q 


26(5) 


i 


j 


2.51 


309 


357 


B1718-32 


0(5) 


— 


E 


— 


3 


T 


2.35 


— 


— 


J1730-2304 


19.9(6) 




E 




3 


T 


0.53 






RI 732—07 

-LJ 1 1 Ji Vf 1 


— 2 4f1 71 


28(3) 


Q 


29(3) 


\ 


x 


2.26 


311 


283 


B1735— 32 


—3(131 




E 




3 


T 


1.20 






RI 7374- 1 3 

i-J \^ -LiJ 


— 22(3) 


— 20(3) 




29(3) 


\ 


x 


1.48 


203 


452 


J1744-1134 


18.72(6) 


-9.5(4) 


c 


20.99(19) 


25 


T 


0.36 


36 


ind. 


B1742-30 


11(4) 


— 


E 


— 


3 


T 


1.91 


— 


— 


B1745— 12 


3(61 




E 




3 


T 


2.45 






B1749— 28 


-4(61 


-5(51 


c 


6(61 


10 


I 


1.23 


35 


9 


RI 7^4—24 


-17(51 




E 




3 


T 


4.40 






B1756-22 


-22(5) 


— 


E 


— 


3 


T 


3.57 


— 


— 


B1757-24 


2(9) 


-3(9) 


C 


4(9) 


26 


I 


5.22 


99 


-12 


RI 802—07 

-LJ lOUii U 1 


-5(61 




E 




3 


T 


3.10 






11 804—271 7 


2(4) 




E 




3 


T 


0.78 






RI 802-1-03 


— 1 2(61 




E 




3 


T 


2.84 






B1804-08 


-5(4) 


1(4) 


C 


5(4) 


10 


I 


2.73 


65 


21 


B1806-21 


-5(7) 


— 


E 


— 


3 


T 


5.23 


— 


— 


B1813— 17 


-9(111 




E 




3 


T 


6.33 






R181 5—14 


— 12(Ql 




E 




3 


T 


7.15 






RI 81 8—04 


— 10(31 


1 0(91 


E 


14(71 


3 


T 


1.94 


129 


34 


B1821+05 


5(11) 


-2(4) 


C 


5(11) 


4 


I 


1.83 


43 


28 


B1820-30A 


2(2) 





E 


— 


3 


T 


8.00 


— 


— 


RI 820—31 


1 6(31 




E 




3 


T 


1.29 






B1821— 24 


—0.90(10) 


—4.6(18) 


c 


4 7(181 


31 


T 


5.80 


129 


ind. 


B1822+00 


— 18(10) 




E 




3 


T 


1.91 






B1822-09 


-13(11) 


-9(5) 


C 


16(10) 


10 


I 


0.88 


67 


9 


B1822-14 


12(14) 




E 




3 


T 


5.06 






B1829-08 


-3(4) 




E 




3 


T 


4.85 






B1829-10 


-13(6) 




E 




3 


T 


5.81 






B1830-08 


-30(3) 


15(11) 


E 


34(6) 


3 


T 


4.66 


751 


163 


B1831-04 


-1(8) 




E 




3 


T 


2.16 






B1834-04 


-2(5) 


12(14) 


E 


12(14) 


3 


T 


4.94 


281 


-18 


B1839+56 


-30(4) 


-21(2) 


C 


37(4) 


4 


I 


1.68 


295 


4 


B1838-04 


6(3) 


-10(8) 


E 


12(8) 


3 


T 


5.68 


323 


-29 


B1842+14 


-9(10) 


45(6) 


C 


46(7) 


4 


I 


2.19 


478 


9 



G. Hobbs et al. 



PSR 


Wong 


Wat 


Coord. 


Wot 


Ref. 


Obs. 


Dist. 


yCL 


AV T 




(mas yr _1 ) 


(mas yr -1 ) 




(mas yr _1 ) 




Type 


(kpc) 


(kms -1 ) 


(kms -1 ) 


B1841-04 


-6(10) 


— 


E 


— 


3 


T 


2.92 


— 


— 


B1841-05 


7(5) 


— 


E 


— 


3 


T 


6.51 


— 


— 


B1844— 04 


-3(6) 




E 




3 


T 


3.26 






Bl 845—01 


— 1(6) 








3 


T 


3.99 






B1846— 06 


17(H) 




E 




3 


T 


3.41 






B1848+13 


-16(6) 


-2(10) 


E 


16(7) 


3 


T 


2.69 


204 


34 


J1852-2610 


-15(11) 


— 


E 


— 


3 


T 


1.75 


— 


— 


Rl 855+02 


-11(71 




B 




3 


T 


7.98 






B1855+09 


-2.94(4) 


-5.41(6) 


c 


6.16(6) 


27 


T 


0.91 


27 


ind. 


B1857— 26 


— 19.9(3) 


—47.3(9) 


c 


51.3(9) 


28 


I 


2.00 


486 




B1859+01 


6(10) 


— 


E 


— 


3 


T 


2.79 


— 


— 


B1900+05 


-8(8) 




E 




3 


T 


4.66 






Bl 900+06 


-4(6) 


— 7(13) 


B 


8(12) 


3 


T 


8.44 


320 


45 


B1900— 06 


-2(8) 




B 




3 


T 


5.37 






B1902— 01 


1(51 


9(14) 


E 


9(14) 


3 


T 


6.01 


256 


38 


B1904+06 


-6(6) 


-1(11) 


E 


6(7) 


3 


T 


8.31 


236 


26 


B1905+39 


11(4) 


11.0(10) 


C 


16(3) 


4 


I 


2.24 


170 


-36 


B1907+00 


-2C61 




B 




3 


T 


3.57 






B1907+02 


1(10) 




E 




3 


T 


4.94 






B1907+10 


-5(A) 


8(8) 


E 


9(8) 


3 


T 


4.18 


178 


g 


J1909-3744 


-9.6(2) 


-35.6(7) 


C 


36.9(7) 


29 


T 


0.46 


80 


86 


J1911-1114 


-6(4) 


-23(13) 


C 


24(13) 


12 


T 


1.22 


139 


42 


B1911+13 


-8(51 


-ill's) 


E 


14(8) 


3 


T 


5.12 


340 


—3 


B1911— 04 


7(13) 


-5(91 


c 


9(12) 


4 


I 


2.79 


119 


18 


B1913+10 


1(31 


-7(6) 


E 


7(6) 


3 


T 


6.27 


208 


— 32 


B1913+16 


-2.56(6) 


0.49(7) 


C 


2.61(7) 


30 


T 


5.90 


73 


15 


B1914+09 


-9(4) 


-1(7) 


E 


9(5) 


3 


T 


2.94 


125 


-3 


B1917+00 




-1(101 


c 




4 


I 


3.06 






B1919+21 




40(10) 


c 




10 


I 


1.09 






B1924+16 


3(81 


— 16(13) 


B 


16(13) 


3 


T 


5.83 


442 


145 


B1929+10 


94.03(14) 


43.4(3) 


C 


103.56(18) 


8 


I 


0.36 


177 


ind. 


B1933+16 


-1(3) 


-13(3) 


c 


13(3) 


10 


I 


5.61 


346 


143 


B1935+25 


-14(31 


-11(4) 


E 


18(4) 


3 


T 


3.25 


277 


-42 


Bl 937+21 


— 0.46(2) 


— 66(2) 


c 


80(2) 


32 


T 


3.60 


14 




B1937— 26 


12(2) 




c 




33 


I 


1.66 


126 


234 


B1944+17 


1(5) 


-9(4) 


c 


9(5) 


6 


1 


1.40 


60 


-23 


B1943-29 


19(9) 




c 




1 


I 


1.54 






Bl 946+35 


— 12.6(6) 


— 7(6) 


c 


12.6(6) 


41 


1 


5.80 


346 


124 


B1951+32 


— 24(4) 


-8(3) 


c 


25 (4) 


34 


1 


2.50 


296 




B1952+29 




—36(10) 


c 




6 


I 


0.70 






B1953+29 


-1.0(3) 


-3.7(3) 


c 


3.8(3) 


19 


T 


4.64 


84 


14 


B1953+50 


-23(5) 


54(5) 


c 


59(5) 


4 


I 


2.24 


626 


-123 


B1957+20 


-16.0(5) 


-25.8(6) 


c 


30.4(6) 


35 


T 


2.49 


359 


-138 


B2000+40 


-13(7) 


-3(7) 


E 


13(7) 


3 


T 


5.90 


364 


148 


B2002+31 


-5(8) 


-10(11) 


E 


11(11) 


3 


T 


7.50 


391 


75 


B2011+38 


-32.1(17) 


-25(3) 


C 


41(2) 


1 


I 


8.44 


1624 


891 


B2016+28 


-2.6(2) 


-6.2(4) 


c 


6.7(4) 


7 


I 


0.97 


31 


ind. 


J2019+2425 


-9.41(12) 


-20.60(15) 


c 


22.65(15) 


36 


T 


1.49 


160 


-62 


B2020+28 


-4.4(5) 


-23.6(3) 


c 


24.0(4) 


7 


I 


2.70 


307 


ind. 


B2021+51 


-5.23(17) 


11.5(3) 


c 


12.6(3) 


7 


I 


2.00 


119 


ind. 



A statistical study of 233 pulsar proper motions 7 



PSR 


Wong 


Wat 


Coord. 


Wot 


Rcf. 


Obs. 


Dist. 


T r/" 1 r 
v T 


AV T 




(mas yr _1 ) 


(mas yr -1 ) 




(mas yr _1 ) 




Type 


(kpc) 


(kms -1 ) 


(kms^ 1 ) 


B2022+50 


17(3) 


14(3) 


E 


22(3) 


3 


T 


2.34 


244 


-56 


B2044+15 


-13(6) 


3(4) 


C 


13(6) 


4 


I 


2.42 


149 


9 


R2043— 04 




— 7(81 


Q 




4 


I 


1.75 






B2045— 16 


117C51 


-5(5) 


c 


117(51 


10 


I 


0.56 


311 


44 


1205 1 —0827 


5.3(10) 


0(3) 


c 


5.3(11) 


37 


T 


1.04 


26 


6 


B2053+21 


-5(6) 


-1(9) 


E 


5(7) 


3 


T 


2.43 


58 


-7 


B2053+36 


-3(7) 


3(3) 


C 


4(6) 


4 


1 


4.62 


88 


18 


B2106+44 


3.5(13) 


1 4,(14) 


c 


3.8(14) 


1 


I 


4.96 


89 


g 


B2110+27 


-23(2) 


— 54(3) 


c 


59(3) 


4 


I 


2.03 


568 


— 179 


B2111+46 


7(12) 


— 11(13) 


E 


13(13) 


3 


T 


4.53 


279 


28 


B2113+14 




-11(5) 


c 


— 


4 


I 


4.17 


— 


— 


J2124-3358 


-14.0(10) 


-47.0(10) 


c 


49.0(10) 


12 


T 


0.27 


63 


-5 


J2129— 5721 


7(2) 


-4C31 


c 


8(3) 


12 


T 


1.36 


52 


45 


J2145— 0750 


— 12.37(6) 


—6.8(8) 


c 


14 1(41 


40 


T 


0.50 


33 


ind. 


B2148+63 


14C31 


10(41 


c 


17C41 


4 


I 


5.51 


444 


656 


B2148+52 


7(3) 


-5(3) 


E 


9(3) 


3 


T 


4.62 


197 


45 


B2154+40 


17.8(8) 


2.8(10) 


c 


18.0(9) 


1 


I 


3.76 


321 


155 


B2217+47 


— 12C81 


-30(6) 


c 


32(7) 


6 


1 


2.22 


337 


35 


B2224+65 


144C31 


112(3) 


c 


182(3) 


4 


1 


1.86 


1605 


121 


J2229+2643 


1(4) 


— 17(4) 


c 


17(41 


19 


T 


1.45 


117 


_2 


B2227+61 


8(10) 


6(11) 


E 


10(11) 


3 


T 


4.04 


192 


79 


J2235+1506 


15(4) 


10(8) 


C 


18(6) 


38 


T 


1.07 


91 


7 


B2255+58 


11(A) 


-14(5) 


E 


18(51 


3 


T 


4.51 


385 


161 


B2303+30 


2(2) 


— 20(2) 


c 


20(3) 


1 


I 


3.66 


347 


25 


B2306+55 


-15(8) 




c 




4 


1 


2.16 






B2310+42 


21.4(9) 


-5.3(12) 


E 


22.0(10) 


3 


T 


1.25 


130 


-31 


J2317+1439 


-1.7(15) 


7(4) 


C 


8(4) 


38 


T 


0.83 


31 


40 


B2315+21 


2(9) 




E 




3 


T 


0.95 






J2322+2057 


-17(2) 


-18(3) 


C 


25(3) 


39 


T 


0.80 


95 


-2 


B2324+60 


-18(5) 


6(5) 


E 


19(5) 


3 


T 


4.86 


438 


-4 


B2327-20 


74.7(19) 


5(3) 


C 


75(2) 


1 


I 


0.39 


138 


36 


B2334+61 


-7(12) 


-8(13) 


E 


11(13) 


3 


T 


3.15 


164 


-35 


B2351+61 


22(3) 


6(2) 


C 


23(3) 


4 


I 


3.43 


374 


-13 



References: 1. Briskcn ct al. (2003a), 2. Freire et al. (2003), 3. Hobbs et al. (2004), 4. Harrison, Lync & Anderson (1993), 5. Siegman, 
Manchester & Durdin (1993), 6. Lync, Anderson & Salter (1982), 7. Brisken et al. (2002), 8. Chatterjee et al. (2004b), 9. van Stratcn ct 
al. (2001), 10. Fomalont ct al. (1997), 11. Wyckoff & Murray (1977) 12. Toscano ct al. (1999b), 13. Splaver et al. (2002), 14. Caravco ct 
al. (1996), 15. Briskcn ct al. (2003b) 16. Dodson et al. (2003), 17. Chatterjee et al. (2001), 18. Langc ct al. (2001), 19. Wolszczan ct al. 
(2000), 20. Bailes et al. (1990a), 21. Konacki, Wolszczan & Stairs (2003), 22. Thorsett et al. (1999), 23. Lewandowski ct al. (2004), 24. 
Camilo, Foster & Wolszczan (1994), 25. Toscano ct al. (1999a), 26. Thorsett, Birsken & Goss (2002), 27. Kaspi, Taylor & Ryba (1994), 

28. Romalont ct al. (1999), 29. Jacoby ct al. (2003), 30. Weisberg & Taylor (2003) 31. Cognard & Kestradc (1997), 32. Cognard et al. 

(1995), 33. McGary ct al. (2001), 34. Migliazzo et al. (2002), 35. Arzoumanian, Fruchter & Taylor (1994), 36. Nice, Splaver & Stairs 
(2001), 37. Doroshcnko ct al. (2001), 38. Camilo, Nice & Taylor (1996), 39. Nice & Taylor (1995), 40. Lohmer et al. (2004). 



The sample contains a wide variety of pulsars including 
'normal', recycled (defined here as pulsars with spin-periods, 
P < 0.1s and spin-down rates, P < 1CF 17 ), 'young' (those 
pulsars with characteristic ages t c < 3Myr) and pulsars 
associated with supernova remnants or globular clusters. In 
our sample, 58% of the proper motions were determined 
using timing methods, 41% using interferometers and the 
remaining 1% using optical telescopes. In Figure 1 we plot 
the positions and proper motions on the sky. The apparent 
motions with respect to the Galactic plane seen in this figure 
are discussed in detail in Section 5. 

Our sample includes 12 pulsars that are situated within 
globular clusters and for which a 2-D speed can be de- 



termined. The mean 2-D speed for this subsample is 
149(8) km s _1 . However, this includes the proper motion 
of the globular clusters as well as the pulsars themselves. 
For 47 Tucanae, removing the Hipparcos measurement of 
the cluster proper motion of /i Q = 5.3(6) mas yr _1 and 
/if = —3.3(6) mas yr _1 (Odenkirchen et al. 1997) from the 
measured proper motions indicates that the intrinsic mean 
2-D speed of the 11 pulsars with measured proper motions 
in the cluster is 25(5) kms -1 . That the mean velocity of re- 
cycled pulsars in globular clusters is significantly lower than 
those moving in the Galactic potential is not surprising given 
the smaller escape velocity of the clusters compared to the 
Galaxy. For 47 Tucanae, for example, the escape velocity is 
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Table 2. Proper motions in Galactic coordinates. In column order, the table lists pulsar name, Galactic coordinates, proper motion in 
Galactic coordinates, the height above the Galactic plane, the velocity in this direction, Vz = fibDcosb, the base-10 logarithms of the 
characteristic age and surface magnetic field strength. 
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~ 58kms _1 (Webbink 1985). As the measured proper mo- 
tion of a pulsar in a globular cluster is dominated by the 
cluster motion we will leave such pulsars out of the remain- 
ing discussion and analysis of this paper. 

In Table 2 we tabulate for each pulsar: Galactic po- 
sition, proper motion in Galactic latitude and longitude, 
height above the Galactic plane, speed perpendicular to the 
Galactic plane, characteristic age, r c = P/2P and surface 
magnetic field strength, B 3 = 3.2 x 10 19 (PP) 1/2 G. The ef- 
fects of stellar motion and Galactic rotation have been re- 
moved from the proper motions given in this table. To do 
this, we assumed a flat Galactic rotation curve with galacto- 
centric distance of the Sun, Rq = 8.5 kpc and the Galactic 
rotation of the Sun, V@ = 225 km s" 1 . These are identical 
to the values used in the analysis by Harrison, Lyne & An- 
derson (1993) and are similar to the current IAU standard 
of Rq = 8.5 kpc and V Q = 220 km s -1 (Kerr & Lynden- 
Bell 1986). Most suggested rotation curves of our Galaxy, 
e.g. Oiling & Merrifield (1998), only deviate significantly 
from a flat rotation curves for distances from the Galactic 
centre < 3 kpc. Since none of the pulsars in our sample are 
within this region of the Galaxy, our assumption of a flat 
rotation curve does not significantly affect our results. 

Similarly, we choose the same values for the peculiar 
motion of the Sun as Harrison, Lyne & Anderson (1993) 
of 15.6 km s" 1 in the direction of / = 48.8° and b = 26.3°. 
This motion is close to the value of 13.4 km s" 1 obtained by 
Dehnen & Binney (1997) using Hipparcos data. Software to 
carry out these conversions between equatorial (or ecliptic) 
and Galactic proper motions is available as part of the ATNF 
pulsar catalogue software package (Manchester et al. 2005) . 
Since the uncertainties in pulsar distances discussed in the 
next section are, on average, much larger than the effects of 
assuming a specific rotation curve and solar motion values, 
the resulting 1-D and 2-D speed distributions are unaffected 
by this correction. In the following, we therefore choose not 
to correct the samples for these small effects. 



3 PULSAR DISTANCES 

Model-independent distances to 22 pulsars within our sam- 
ple are available from measurements of annual parallax. 
These distances are known to within 10% accuracy and are 
listed in Table 3. Limits on the distances for a further 13 pul- 
sars have been determined using HI absorption or by asso- 
ciations with globular clusters or supernova remnants (such 
limits are reviewed by Frail & Weisberg 1990). The distances 
to the remaining pulsars are estimated from their disper- 
sion measures and a Galactic electron density model. The 
TC93 model attempted to predict distances accurate to 25% 
or better on average. The CL02 model purports to provide 
distances with a similar accuracy (Cordes & Lazio 2002b 
report that typical distance uncertainties are slightly less 
than 20%). Extensive comparisons between the two models 
have been made (see, for example, Cordes & Lazio 2002b 
and Kramer et al. 2003). Here, we emphasise that the two 
models are not entirely consistent and list major discrepan- 
cies for the pulsars considered in this paper in Table 4. With 
the CL02 model, only a lower limit of 50 kpc is provided for 
the distance to PSR B1620-09 (the TC93 model places this 
pulsar at 3.9 kpc). The CL02 distance is clearly incorrect 
and hence, we do not include this pulsar in the remaining 
analyses presented in this paper. 

As the distances predicted by the CL02 electron den- 
sity model are, on average, smaller than the TC93 model, 
the derived 2-D speeds are also lower (Figure 2a). For our 
sample, the mean distance has decreased by ~2% (the me- 
dian decrease is ~10%). According to the TC93 model, four 
pulsars have 2-D speeds V^f c > 1000 kms -1 , whereas, us- 
ing the CL02 model, only two pulsars (PSRs B2011+38 
and B2224+65) have such large 2-D speeds. The appar- 
ent 2-D speeds of the remaining two, PSRs B0523+11 and 
B2148+63, have decreased to V$ L = 446 and 449 kms" 1 
respectively. 

In passing, we note that the CL02 model places pul- 
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sars at distances almost halfway between the Lyne, Manch- 
ester & Taylor (1985) model (that was commonly used prior 
to the TC93 model) and the TC93 model (Figure 2b). For 
the remainder of this paper we use the CL02 model (un- 
less independent distance estimates exist), but caution that, 
even with this most recent distance model, the distances to 
some pulsars, and hence their velocities, may be significantly 
under- or over-estimated. We therefore emphasise the im- 
portance of continuing to measure pulsar distances to refine 
future electron density models. 



4 ONE AND TWO-DIMENSIONAL SPEEDS 

For all pulsars in our sample at least one component of their 
space velocity is determined. We make use of these by defin- 
ing a pulsar's one-dimensional speed Vi = fiD, where fi is the 
proper motion in either the longitudinal or latitudinal direc- 
tion and D is the distance obtained using the CL02 model. 
For 156 pulsars we can obtain their 2-D speeds Vi = jitotD, 
where /itot is the total proper motion. In what follows, we 
investigate aspects of the 1-D and 2-D speeds obtained for 
the sample as a whole and in various subsets of the sample. 

4.1 Statistical overview 

A summary of the 1-D and 2-D speeds for various subsets of 
the complete sample is given in Table 5. Histograms of the 
data are presented in Figures 3 and 4. Table 5 contains a 
description of each subset, the number of pulsars within this 
set, the number of measured velocity components (for the 
1-D speeds, this is less than twice the number of pulsars as 
some pulsars have only a proper motion measurement in one 
coordinate) and average distances and characteristic ages. 

The mean 1-D speed is 133(8) kms -1 and the fastest 
moving pulsar in our sample, PSR B2011+38, has V\ = 



o 



Z o 
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0.1 1 10 100 1000 10 

1-D speed (km s~ 1 ): t q > 3Myr 




0.1 1 10 100 1000 10 



1 — D speed (km s ): r c < 3Myr 

Figure 3. One dimensional speeds for pulsars in our sample that: 
a) have characteristic ages greater than 3 Myr; (b) have charac- 
teristic ages less than 3 Myr. 
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Table 3. Pulsars with parallax measurements. This table con- 
tains each pulsar's name, parallax, distance as measured from 
the parallax, a reference for the measurement and the pulsar's 
2-D speed assuming the proper motion given in Table 1 and the 
distance determination. 



pen 
rill 


Parallax 


uist. 


net. 






(mas) 


(kpc) 




^Km S ) 


B0329+54 


0.94(11) 


1.06(12) 


1 


98(11) 


J0437-4715 


6.98(19) 


0.143(4) 


2 


94(3) 


J0633+1746 


6.4(18) 


0.16(5) 


3 


125(40) 


B0656+14 


3.47(36) 


0.29(3) 


4 


59(7) 


B0809+74 


2.31(4) 


0.433(7) 


1 


103(2) 


B0823+26 


2.8(6) 


0.36(8) 


5 


186(42) 


B0833-45 


3.4(2) 


0.29(2) 


6 


80(5) 


B0919+06 


0.83(13) 


1.20(19) 


7 


503(80) 


B0950+08 


3.82(7) 


0.262(5) 


1 


36(1) 


B1133+16 


2.80(16) 


0.36(2) 


1 


641(36) 


B1237+25 


1.16(8) 


0.86(6) 


1 


481(33) 


B1257+12 


1.3(4) 


0.8(2) 


8 


344(90) 


B1451-68 


2.2(3) 


0.45(6) 


9 


88(12) 


J1713+0747 


0.9(3) 


1.1(4) 


10 


34(13) 


J1744-1134 


2.8(3) 


0.36(4) 


11 


36(4) 


B1855+09 


1.1(3) 


0.9(2) 


12 


27(6) 


B1857-26 


< 1.1 


> 0.9 


13 


>219 


B1929+10 


3.02(9) 


0.331(10) 


1 


177(5) 


B2016+28 


1.03(10) 


0.97(9) 


1 


31(3) 


B2020+28 


0.37(12) 


2.7(9) 


1 


310(100) 


B2021+51 


0.50(7) 


2.0(3) 


1 


120(18) 


J2145-0750 


2.0(6) 


0.5(2) 


14 


33(9) 



References: 1. Briskcn ct al. (2002), 2. van Straten et al. (2001), 
3. Caraveo ct al. (1996), 4. Briskcn ct al. (2003b), 5. Gwinn ct 
al. (1986) 6. Dodson et al. (2003) 7. Chattcrjcc ct al. (2001) 8. 
Wolszczan et al. (2000) 9. Bailcs ct al. (1990b) 10. Camilo, Foster 
& Wolszczan (1994) 11. Toscano ct al. (1999a) 12. Kaspi, Taylor 
& Ryba (1994) 13. Fomalont ct al. (1999) 14. Lohmer et al. (2004) 

Table 4. The ten largest distance changes in our sample between 
the TC93 and CL02 models. 



PSR 


d tc 


D CL 




(kpc) 


(kpc) 


B0523+11 


>7.7 


3.1 


B0559-05 


>7.5 


3.9 


B0736-40 


>11.0 


2.6 


B0835-41 


4.2 


1.1 


B1552-31 


>6.1 


2.3 


B1737+13 


>4.8 


1.5 


B1802-07 


11.6 


5.0 


B1900-06 


8.8 


5.4 


B2011+38 


13.1 


8.4 


B2148+63 


>13.7 


5.5 



1284 km s 1 in right ascension and 996 km s 1 in declination. 
The smallest upper limit is V\ < 2kms _1 in right ascension 
for PSR J1640+2224. The 2-D speeds range from 12kms _1 
for PSR B1706-16 to 1625 kms" 1 for PSR B2011+38. 
The mean 2-D speed is 211(18) kms -1 , or 269(25)kms~ 1 
with the TC93 model. This latter value is consistent with 
300(30) kms -1 obtained by Lyne & Lorimer (1994) who also 
used the same distance model. 

4.2 Isotropy of the velocity vector 

The observed 1-D and 2-D speeds we measure are, of course, 
the projection of the pulsar's 3-D space velocity. For an 
isotropically distributed pulsar velocity vector of magnitude 
V3, it is straightforward to show that the mean 1-D and 2-D 
speeds are given by V3/2 and 71-V3/4 respectively. A simple 
test of isotropy, then, is to compare the ratios of the mean 2- 
D and 1-D speeds which should be (V 2 )/(Vi) = tt/2 ~ 1.57. 
For the sample as a whole, we find (V2) / (Vi) = 1.6(2) in ex- 
cellent agreement with an isotropic velocity vector. Indeed, 
all the subsamples in Table 5 give similar results. This simple 
but important result will be used to derive the 3-D velocity 
distribution in Section 6. 

4.3 High-velocity pulsars 

The highest 2-D speeds in our sample are 1624, 1608 and 
740 kms" 1 for PSRs B2011+38, B2224+65 and B1830-08 
respectively 4 . As there is a large gap between the second 
and third fastest objects, we must consider whether the 
speeds greater than 1000 km s^" 1 are real or due to inaccu- 
racies either in the distance model or in the proper mo- 
tion measurements. If the true 2-D speed of PSR B2011+38 
were lOOOkms -1 then the observed proper motion would 
be 25masyr~ 1 for its estimated distance of 8 kpc. This 
required proper motion is significantly lower than recent 
VLA measurements by Brisken et al. (2003) who obtain 
iitot = 40(2) masyr -1 and from Paper I where we measured 
Utot = 51(10) masyr -1 . We can therefore rule out the proper 
motion measurements being in error. However, as there are 
no independent distance estimates to either PSR B2011+38 
or B2224+65 it is possible that the reported distances may 
be overestimated (the reported distance to PSR B2011+38 
of 8.4 kpc is the largest in our sample). 

If these 2-D speeds are correct, then it means that the 
fastest pulsars are still part of a long tail in the velocity dis- 
tribution with the gap between them and the slower-moving 
pulsars being caused by small number statistics (i.e. pulsars 
at the high end of the distribution are rare and so the gap 
is only an artifact) and do not make up a separate popu- 
lation of high-velocity objects. To investigate whether this 
is possible, we undertook a simple simulation where we ap- 
proximated Figure 2 by a single 2-D speed distribution. In 
three out of fifteen Monte Carlo realisations of the sample, 
we found a gap above and below 1000 km s -1 . We there- 
fore conclude that the gap observed in Figure 2 is entirely 



4 These are the observed 2-D speeds of the pulsars. After remov- 
ing the effects of Galactic rotation and the peculiar motion of the 
Sun, the 2-D speeds in the pulsar's local standard of rest become 
1459, 1640 and 730 km s" 1 respectively. 
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Table 5. Parameters for different samples of pulsars. The top half of the table provides values for 1-D speeds, the bottom half contains 
2-D speeds. Np Sr and Np Sr indicate the number of pulsars for which 1-D and 2-D speeds have been measured. The Np Sr pulsars provide 
Ng 1-D velocity components. The table provides mean and median speeds (V, V med ), distances (D, D med } an d characteristic ages (f c , 
T c " le<! ) and an estimate of the mean space velocity V3 obtained from the 1-D and 2-D speeds; see text. 





All 


Normal 


Recycled 


SNR 


r c < 3 Myr 


r c > 3 Myr 


D < 500 pc 


Ml 


217 


178 


39 


g 


73 


105 


19 


N S 


372 


299 


73 


15 


119 


180 


36 


Vi (kins- 1 ) 


133(8) 


152(10) 


54(6) 


150(42) 


192(20) 


126(10) 


77(17) 


yrned ( kms -l) 


79 


100 


38 


95 


123 


85 


54 


D (kpc) 


2.36(12) 


2.58(14) 


1.32(19) 


2.4(6) 


3.3(3) 


2.06(13) 


0.35(2) 


D med ( kp( .) 


1.91 


2.14 


0.95 


2.25 


2.79 


1.83 


0.36 


f c (Myr) 


1376 


70.7 


7333 


0.1 


1.2 


118.9 


1719 


T med (Myr) 


7.4 


4.3 


5610 


0.07 


1.1 


14.1 


37 


V 3 = 2V 1 (kms- 1 ) 


266(16) 


304(20) 


108(12) 


300(82) 


384(40) 


252(20) 


154(34) 


N 2 

psr 


156 


121 


35 


7 


46 


75 


17 


V 2 (kms" 1 ) 


211(18) 


246(22) 


87(13) 


227(85) 


307(47) 


209(19) 


128(33) 


V 2 med (kms- 1 ) 


139 


178 


73 


142 


240 


152 


94 


D 2 (kpc) 


2.25(15) 


2.50(17) 


1.4(2) 


2.7(7) 


3.0(3) 


2.18(17) 


0.35(2) 


D med ( kpc ) 


1.70 


2.07 


1.04 


2.5 


2.10 


1.99 


0.36 


fc,2 (Myr) 


1732 


59.8 


7511 


0.06 


1.1 


95.7 


1565 


T med ( Myr ) 


8.7 


4.82 


5210 


0.04 


0.9 


10.2 


18 


V 3 = (4/n)V 2 


269(23) 


313(28) 


111(17) 


289(108) 


391(60) 


266(24) 


163(42) 
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Figure 4. Two-dimensional speeds for: (a) all the pulsars in our sample; (b) pulsars with characteristic ages less than 3 Myr; (c) recycled 
pulsars; (d) normal (not recycled) pulsars with characteristic ages greater than 3 Myr. The mean 2-D speeds for each sample are 211(18), 
307(47), 87(13) and 209(19) kms" 1 respectively. Panels b, c and d arc drawn to the same scale. 
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consistent as being the result of small-number statistics of 
sampling a continuous velocity distribution with a large tail. 

To confirm these high 2-D speeds will require indepen- 
dent distance measurements. However, the observed bow 
shock for PSR B2224+65 does suggest a high velocity for 
this pulsar which is most likely greater than 800 km s~ (see 
Cordes, Romani & Lundgren 1993 and Chatterjee & Cordes 
2004a). 

4.4 Young and old pulsars 

After selecting only those pulsars with characteristic ages 
less than 3 Myr (Figure 4b), we obtain a mean 2-D speed of 
307(47) kms -1 , consistent with the value of 345(70) kms -1 
obtained by Lyne & Lorimer (1994) who estimated distances 
using the TC93 model. The mean 2-D speed for young pul- 
sars is higher than the value of 209(19) km s -1 for older, 
non-recycled pulsars; see Figure 4d. However, the veloc- 
ity distributions for the young and older pulsars appear to 
overlap. We note that the sample of 46 young pulsars con- 
tains nine objects with 2-D speeds less than 100 km s -1 and 
only seven with 2-D speeds > 400 kms -1 . The smallest 2-D 
speed, 12kms -1 , is measured for PSR B1706— 16, although 
whether this represents a truly low-velocity population or a 
projection effect is not clear. We investigate this issue fur- 
ther in Section 6. 

4.5 Pulsar — supernova remnant associations 

It is notable that the 1-D and 2-D speeds for the eight pul- 
sars that have been associated with supernova remnants (see 
Tables 5 and 6) are significantly lower than the mean for 
other young pulsars and consistent with the mean for the 
entire pulsar sample that includes both young and old pul- 
sars. In contrast, by inferring 2-D speeds based on the offset 
of pulsars from the centre of the supernova remnant, Frail, 
Goss & Whiteoak (1994) suggested that most pulsars in su- 
pernova remnants have velocities that are significantly larger 
than the mean. More recent proper motion measurements 
show that some of the velocities claimed by Frail et al. are 
too high. For example, PSR B1757— 24 was thought to be in 
G5.4— 1.2, giving rise to a 2-D speed which was predicted to 
be 1800 kms -1 but which has subsequently been shown to 
be ~100kms _1 , so that it is unlikely to be associated with 
the supernova remnant (Gaensler & Frail 2000; Thorsett, 
Brisken & Goss 2002). Other pulsars with 2-D speeds pre- 
dicted by Frail et al. to be greater than 1000 kms -1 arc 
PSRs B1509-58, B1610-50, B1800-21 and B1930+22. Un- 
fortunately these pulsars glitch regularly and it is difficult 
to obtain proper motions using timing techniques. Interfer- 
ometric measurements of these pulsars are therefore highly 
desirable. 

One explanation for the observed low velocities of pul- 
sars associated with supernova remnants is that any high- 
velocity objects would leave the central regions of the rem- 
nant within a short timescale and are therefore hard to de- 
tect in a targeted search of a supernova remnant. The dif- 
ficulty in identifying a genuine association is additionally 
compounded by the fact that the probability of a chance 
association with a pulsar increases with the square of the 
angular offset from the remnant centre. One possible ex- 
ample of a high-velocity pulsar associated with a supernova 



Table 6. Pulsars in our sample that have been associated with 
supernova remnants. Recent work (e.g. Thorsett, Brisken & 
Goss 2002) indicates that PSR B1757— 24 is not associated with 
G5.4-1.2 as was previously thought. For PSR J0538+2817 only 
one component of velocity has been determined. 



PSR 


SNR 


Tc 


V T 






(kyr) 


(kms -1 ) 


B0531+21 


Crab 


1.24 


140 


J0538+2817 


S147 


618 


>140 


B0656+14 


Monogcm Ring 


111 


59 


B0833-45 


Vela 


11.3 


80 


(B1757-24) 


(G5.4-1.2) 


(15.5) 


(89) 


B1830-08 


W41 


41 


740 


B1951+32 


CTB80 


10.7 


300 


B2334+61 


G114.3+0.3 


40.9 


159 



remnant is PSR B1830— 08 (nominally the third fastest ob- 
ject in our sample with Vf L = 740 kms -1 ) and the shell 
supernova remnant W41 (Clifton & Lyne 1986; Gaensler 
& Johnston 1995). Regardless of the distance uncertainties, 
the new proper motion measurement for PSR B1830— 08 
presented in paper I shows that the pulsar is moving away 
from W41 and is consistent with birth within the bound- 
aries of the remnant. Further observations to constrain the 
distance and origin of this young pulsar are required. 

In Table 6 we include an association between 
PSR B2334+61 and G114.3+03 which was first suggested by 
Kulkarni et al. (1993) and Ffirst, Reich & Seiradakis (1993). 
Fiirst et al. found that the ages and distances of both objects 
are comparable and that the pulsar is located close to the 
centre of the remnant. They also noted that, assuming the 41 
kyr characteristic age to be a good indicator of the pulsar's 
true age, the expected proper motion from birth in the cen- 
tre of the supernova remnant to the pulsar's current position 
should be /itot = 12(2) mas yr -1 . Our timing measurement 
of fitot = 11(13) masyr -1 is consistent with the predicted 
value for an association. Clearly, given the large uncertainty 
in utot, we cannot rule out that the pulsar was not born in 
the centre of the supernova remnant. We can state, however, 
that B2334+61 was born somewhere within G114.3+0.3 
(even with a proper motion of fitot = 24 masyr -1 , the pul- 
sar would only have moved ~16arcmin or only 25% of the 
current diameter of the supernova remnant within 41 kyr). 
The birth position of B2334+61 will be determined more ac- 
curately with a few more years of timing data or by making 
interferometric measurements. 



4.6 Recycled pulsars 

In Table 7 we list the 2-D speeds for the recycled pulsars 
in our sample. The mean value of Vr = 87(13) kms -1 is 
entirely consistent with the value of 85(13) kms - reported 
by Toscano et al. (1999) using the TC93 model. This con- 
sistency is due to the known millisecond pulsars being rel- 
atively close-by and so CL02 and TC93 agree much better 
than they do for more distant pulsars. Also, a parallax has 
been measured for many of these pulsars providing indepen- 
dent distance estimates (for the millisecond pulsars, only 
the predicted distances to PSRs J0218+4232, J1045-4509, 
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Table 7. 2-D speeds for the recycled pulsars. For binary systems, 
the orbital period Pi is also given. The orbital period listed for 
PSR B1257+12 is that of the most massive planet; see text. 



PSR 


Vt 


Pi, 




(kms -1 ) 




J0034-0534 


79 


1.59 


J0218+4232 


63 


2.03 


Tfl/1 QV A 71 
JU4o f — 4 t 10 


Q/1 

y4 


^ 7/1 
D. ( 4 


J0613-0200 


59 


1.19 


J0621+1002 


23 


8.31 


J0711-6830 


89 


— 


J0751+1807 


>1+ 


0.26 


J1012+5307 


49 


0.60 


J1022+1001 


>35t 


7.81 


J1024-0719 


150 


— 


J1045-4509 


72 


4.08 


B1257+12 


345 


66.5 


Ti ii r r QQQO 

J1400— oooU 


63 


76.2 


J1518+4904 


28 


8.63 


B1534+12 


113 


0.42 


J1603-7202 


47 


6.31 


J1640+2224 


73 


176 


J 104o — 1ZZ4 




1/17 
14 ( 


J1709+2313 


68 


22.7 


J1713+0747 


34 


67.8 


J1730-2304 


>50t 


— 


J1744-1134 


36 




T1 CO/1 T71 *7 


> ( ' 


n.i 


B1855+09 


27 


12.3 


J1909-3744 


80 


1.53 


J1911-1114 


139 


2.72 


B1913+16 


73 


0.32 


B1937+21 


14 




B1953+29 


84 


117 


B1957+20 


359 


0.38 


J2019+2425 


160 


76.5 


J2051-0827 


26 


0.10 


J2124-3358 


63 




J2129-5721 


52 


6.63 


J2145-0750 


33 


6.84 


J2229+2642 


117 


93.0 


J2235+1506 


91 




J2317+1439 


31 


2.46 


J2322+2057 


95 





f For these pulsars, only the proper motion coordinate in ecliptic 
longitude has been measured. The given value is therefore a lower 
limit to its 2-D speed. 

J1643-1224, J2129-5721 and J2317+1439 have changed by 
more than 1 kpc from the TC93 to the CL02 model). 

Toscano et al. (1999) also noted that isolated millisec- 
ond pulsars have lower velocities than binary millisecond 
pulsars (see also Johnston et al. 1998). Our sample con- 
tains seven solitary and 28 binary millisecond pulsars with 
measured 2-D speeds. The size of this sample more than 
doubles that of Toscano et al. (1999). The mean 2-D speeds 
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Figure 5. 2-D velocity versus orbital period for binary millisec- 
ond pulsars. 



we obtain for the solitary and binary millisecond pulsars are 
77(16) and 89(15) kms~ respectively which are not signifi- 
cantly different. 

Tauris & Bailes (1996) predicted that the velocities of 
binary millisecond pulsars should increase with decreasing 
orbital period. No such correlation was found in their data 
which they attributed to selection effects. Our data, shown 
in Figure 5, also show no correlation. Tauris & Bailes (1996) 
predict a lower limit for the velocity of such pulsars of 
« 25kms _1 and an upper limit of ~270kms~ 1 . Two pul- 
sars in our sample, PSRs B1257+12 and B1957+20, with 2- 
D speeds of 344 and 359 km s _1 respectively, are clearly well 
above this upper limit. The first, PSR B1257+12 is known to 
be part of a planetary system and may have undergone a dif- 
ferent evolution than other binary millisecond pulsars. The 
second, the 'black-widow' pulsar PSR B1957+20, is thought 
to be slowly ablating its companion (Fruchter, Stinebring 
& Taylor 1988) and may, therefore, provide a link between 
binary and isolated millisecond pulsars. This seems rather 
unlikely since we do not see solitary millisecond pulsars with 
such high velocities. 

We may expect double neutron star systems to have 
even lower velocities than millisecond pulsars as a conse- 
quence of having survived two supernova explosions. Un- 
fortunately, we currently only have three double neutron 
star systems with proper motion determinations. Our proper 
motion measurement for PSR J1518+4904 from Paper I of 
25(7) kms -1 agrees well with 28(2) kms -1 found by by Nice, 
Taylor & Sayer (1999) (only the uncertainties in the total 
proper motion were considered when determining the errors 
on these values). Earlier measurements for PSR B1534+12 
and B1913+16 imply respective 2-D speeds of 113.0(1) and 
73(2)kms _1 respectively. More recently, from scintillation 
measurements of the double pulsar J0737— 3039, Coles et 
al. (2004) derive a 2-D speed of ~66kms _1 . A proper mo- 
tion measurement of this system should be available soon 
through both interferometric and timing measurements. 
Given the small-number statistics, we simply note that these 
2-D speeds are entirely consistent with the distribution for 
the whole sample of recycled pulsars. 
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5 MOTION IN THE GALACTIC PLANE 

Standard models of pulsar birth suggest that pulsars are 
born in Type II supernovae (see, for example, Gunn & 
Ostriker 1970). The progenitor massive stars lie along the 
Galactic plane with a scale height of ~63pc (Sun & Han 
2004). Any velocity kick imparted to the pulsar at its birth 
will make it move away from the Galactic plane. Hence, the 
motion of young pulsars that appear to be moving toward 
the plane can only be explained by projection effects or by 
the pulsar being born well above the plane and out of the 
expected progenitor region. For instance, the unknown ra- 
dial motion of a pulsar may make it appear to be moving 
towards the plane when it is, in fact, moving away (see, for 
example, Helfand & Tademaru 1977). Restricting our study 
to pulsars with characteristic ages less than 3 Myr guaran- 
tees that all bound objects are within the first quarter-cycle 
of their oscillation perpendicular to the Galactic plane (see, 
for example, Cordes & Chernoff 1998). 

Harrison, Lyne & Anderson (1993) found that the po- 
sitions and proper motions of PSRs B0450-18, B0523+11, 
B0540+23 and B1718-02 could only be explained by the 
pulsars being born above the Galactic plane and currently 
moving moving toward the plane. However, Cordes & Cher- 
noff (1998) analysed the proper motions of 49 pulsars with 
characteristic ages less than 10 Myr and found that all were 
consistent with the earlier conclusion of Helfand & Tade- 
maru (1977) that all pulsars are born in the plane and sub- 
sequently move away. The main difference between the anal- 
yses was that Harrison, Lyne & Anderson (1993) defined 
a much smaller progenitor region than Cordes & Chernoff 
(1998). For instance, even though PSR B0540+23 is likely to 
be moving towards the Galactic plane, due to its low height 
above the Galactic plane of only 0.2 kpc, it was deemed, by 
Cordes & Chernoff (1998), to have been born within the 
expected pulsar progenitor region of the Galaxy. 

In our sample, only PSR B0114+58 (and B0540+23; see 
above discussion) has a characteristic age less than 3 Myr, a 
current z-height greater than 100 pc and is seemingly moving 
toward the plane. However, this pulsar lies below the plane 
at a distance of 130 pc which is within current estimates of 
the pulsar progenitor scale height. We therefore agree with 
Cordes & Chernoff (1998) that all the pulsars that seem to 
be moving toward the plane are consistent with having been 
born within the plane. 

In Figure 6 we indicate the angle of velocity vectors with 
respect to the Galactic plane for pulsars with Galactic lati- 
tudes less than 30° . The lengths of the tracks are indicative 
of the magnitude of the pulsar velocities and were calculated 
using the pulsar proper motions in Galactic longitude and 
latitude only; as before, the radial velocity component has 
been ignored. We find no predominant angle for the motion 
of the pulsars with respect to the Galactic plane. 

Earlier we analysed pulsar velocities based on proper 
motion measurements in celestial or ecliptic coordinates. 
The difference between these velocities and those obtained 
after the effects of Galactic rotation and the peculiar mo- 
tion of the Sun have been removed is minimal and is smaller 
than the likely uncertainties due to the pulsar distances. The 
shape of the 2-D speeds after taking into account the effects 
of Galactic rotation is little changed from the earlier version 
(the mean velocities are identical). However, the velocities 
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Figure 6. The motion of pulsars with respect to the Galactic 
plane. Pulsars indicated by a circle have r c < 1 Myr, those with 
t c > 1 Gyr by a '*' symbol. All other pulsars are indicated using 
an 'x'. 



for a few pulsars change significantly. For instance, in the 
most extreme case the velocity of PSR B1933+16 reduces 
from 346 to 208 km s -1 . 



6 THREE-DIMENSIONAL VELOCITIES OF 
YOUNG PULSARS 

With the exception of two binary pulsars (Bell et al. 1995; 
van Kerkwijk et al. 1996) it is generally not possible to make 
a direct measurement of the true 3-D space velocity of a 
pulsar. However, as the observed data are consistent with an 
isotropic velocity vector, it is possible to determine the 3-D 
speed distribution from the observed 1-D and 2-D speeds. 
In this section we describe and discuss the consequences of 
a new technique that deconvolves the expected form of a 
single-magnitude isotropic space velocity from the observed 
1-D and 2-D distributions. 



6.1 The deconvolution technique 

Following Lyne & Lorimer (1994), we derive the birth ve- 
locity distribution using the sample of pulsars with r c < 3 
Myr to minimise any selection biases against detecting older 
high-velocity pulsars. The resulting samples of 119 1-D and 
46 2-D speeds therefore represent two projected views of 
the underlying 3-D distribution which we seek to recover. 
For a single-magnitude isotropic space velocity distribution 
of Vz, the probability density function corresponding to a 
logarithmic bin Vi for 1-D speeds 

P(Vi)dlog(Vi) = §-dlog(Vi) 

V3 



(1) 



(0 < Vi < V 3 ). 
Similarly, for the 2-D case we find 



P(V2)dlog(V 2 ) = 



vi d\og{v 2 ) 
Vz {vi - v?y/* 
(0 < v 2 < v z ). 



(2) 



One way to think about these 1-D and 2-D distributions 
is that they represent the convolution of the underlying 3- 
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Figure 7. The 3-D velocity distributions obtained from the observed 1-D (left) 2-D (right) distributions using the deconvolution technique 
described in the text. The uncertainties on each histogram bin are calculated as y/N where N is the number of pulsars in each bin. The 
histograms arc normalised to unit area. The dotted curve shows the 3-D distribution favoured by Arzoumanian ct al. (2002). The solid 
curve is the best-fitting Maxwcllian distribution to the histogram from the 2-D distribution with a = 265 km s _1 . 



D distribution with the 1-D and 2-D probability density 
functions. We therefore can use a CLEAN-type algorithm 
(Hogbom 1974) to deconvolve these functions from the ob- 
served speed distributions. In detail, we (i) cross-correlate 
the analytic expression with the observed distribution to 
determine the velocity, V3, corresponding to the maximum 
correlation; (ii) subtract a scaled version of the analytic ex- 
pression defined by V3 from the observed distribution with a 
small gain factor; (iii) record a "CLEAN" -component corre- 
sponding to V3 and repeat for a large number of iterations. 
The "CLEAN" -components are subsequently binned to give 
a histogram representing the true 3-D space velocites. 

6.2 The underlying space velocity distribution 

The resulting 3-D distributions obtained from the 1-D and 
the 2-D observed distributions are shown in Figure 7. Both 
forms of the 3-D distribution are in excellent agreement 
with each other and have mean 3-D speeds of 400(40) and 
431(60) kms -1 respectively. The smaller error on the mean 
derived from the 1-D speeds reflects the larger size of this 
sample. It is notable that the long "tail" towards low speeds 
in the observed 1-D histogram (Figure 3b) is not present in 
the 3-D histogram. This reflects the fact that the observed 
"tail" is simply due to projection effects. Our 3-D distribu- 
tions show that very few young pulsars have 3-D speeds less 
than ~60kms -1 . The right panel of Figure 7 shows that an 
excellent fit (reduced \ 2 = 0.6) can be obtained using a sim- 
ple Maxwellian distribution with a 1-D rms a = 265 km s _1 . 

We have tested our deconvolution algorithm by simulat- 
ing 3-D velocity distributions of different forms, randomly 
selecting as many 1-D speeds as in our real sample from 
these 3-D distributions, binning the resulting 1-D speeds as 
in Figure 3 and using our deconvolution method to derive 
the underlying 3-D velocities. In all cases, we have success- 
fully reproduced the original form of the simulated distribu- 
tion. It is of particular interest to determine whether a sin- 
gle component velocity distribution could be obtained using 
our algorithm from a simulated bimodal 3-D distribution. 
We have therefore formed a bimodal 3-D distribution from 
two Gaussians of differing amplitudes centred on 100 and 
630 km s -1 similar to that proposed by Arzoumanian et al. 
(2002). Using the technique described above, we attempted 



to derive a 3-D distribution using our deconvolution method. 
In all cases we reproduced a bimodal distribution. The zero- 
lag cross correlation between the predicted and the simu- 
lated distributions provided a measure of the goodness of 
the reproduced distribution. All the correlation coefficients 
were found to be close to 1 implying that the deconvolution 
techqniue accurately reproduces the bimodal distribution, 
the positions of the components and their relative ampli- 
tudes. Similar results are obtained when simulating a single 
component distribution. 

6.3 Comparison with earlier work 

Although most authors agree that a broad spectrum of ve- 
locities is required to explain the observations, considerable 
disagreement exists on whether the distribution has more 
than one component. Based on a number of available con- 
straints, Fryer, Burrows & Benz (1998) favour a bimodal 
velocity distribution with one peak near Okms -1 and an- 
other above 600 km s -1 . Similarly, from an extensive study 
of a sample of pulsars detected in low-frequency (~ 0.4 GHz) 
radio surveys, Arzoumanian, Chernoff & Cordes (2002) find 
a two component distribution with characteristic velocities 
of 90 and 500kms _1 . Whether the origin of these two com- 
ponents is due to a physical effect of the kick mechanism is 
still a matter for debate. 

The results of our novel deconvolution analysis of a large 
sample of pulsars do not support the idea of a bimodal distri- 
bution of pulsar velocities. This can be seen in the left panel 
of Figure 7 where we compare our results with an appro- 
priately normalised version of the 3-D distribution preferred 
by Arzoumanian et al. (2002) . This figure clearly shows that 
the new distribution is incompatible with the Arzoumanian 
et al. model which predicts an excess of low- velocity pulsars 
and dearth of high- velocity ones. Fryer et al. (1998) argued 
that the pulsar proper motion data available to them at that 
time were not sufficient to constrain distribution of neutron 
star kick velocities. It would be instructive (but beyond the 
scope of this paper) to repeat their calculations using the 
larger sample of proper motions now available. For now, we 
re-iterate that a Maxwellian distribution provides an excel- 
lent fit to the 3-D pulsar velocity distribution. 

An earlier analysis by Hansen & Phinney (1997) on a 
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smaller sample of pulsars also preferred a Maxwellian distri- 
bution for V3. However, their value of a = 190 km s -1 does 
not provide a very good fit (reduced \ 2 — 5.7) to the new 
larger sample presented here. Our results, however, do agree 
well with the mean velocity of 450(90) km s -1 found by Lyne 
& Lorimer (1994). The larger sample of pulsars in this anal- 
ysis has resulted in a reduction of the uncertainties by more 
than a factor of two. The slight reduction in the mean value 
of 400(40) km s -1 compared with Lyne & Lorimer (1994) 
is largely a result of the smaller average distances inferred 
from the CL02 model compared with TC93. 

The different pulsar velocity distributions proposed pre- 
dict different fractions of pulsars that will escape the Galac- 
tic gravitational potential. Assuming the escape velocity of 
our Galaxy to be ~ 430 kms -1 (Leonard & Tremaine 1990), 
we find that 35% of all pulsars will escape the Galactic po- 
tential in the new distribution compared to 50% for the bi- 
modal distribution of Arzoumanian et al. (2002) and only 
16% for the Hansen & Phinney (1997) distribution. 

In summary, our new results on the 3-D pulsar velocity 
distribution more closely match the earlier conclusions of 
Lyne & Lorimer (1994) than the results of Hansen & Phin- 
ney (1997), Fryer, Benz & Burrows (1998) or Arzoumanian 
et al. (2002). While we have tried to account for selection 
effects by placing a 3 Myr cut-off in our sample of young pul- 
sars, a more detailed analysis that fully takes into account 
the various selection effects present in this new sample (and 
ideally including the other constraints discussed by Fryer et 
al. 1988) should now be carried out to further investigate 
the results presented here. 



The implications of these results for 'kick' mechanisms 
may be summarised by stating that the true space veloc- 
ities of young pulsars range from a few 10s to well over 
1000 kms -1 with a mean velocity of 400(40) kms -1 . Ac- 
cording to Lai, Chernoff & Cordes (2001): (1) local con- 
vective instabilities in the collapsed stellar core can account 
for velocities up to ~ 100 kms -1 ; (2) global asymmetric per- 
turbations can create velocities over 1000 kms -1 ; (3) asym- 
metric neutrino emission can provide kick velocities up to 
~ 1000 kms -1 ; (4) the electromagnetic rocket effect can ac- 
celerate pulsars up to similarly high velocities. Our results 
suggest that (1) is unlikely. (4) is testable by studying the 
alignment between the direction of motion of a pulsar and 
its spin-axis (see Deshpande, Ramachandran & Radhakrish- 
nan 1999) although the duration of the kick will also affect 
the observed alignment (Spruit & Phinney 1998). 

More proper motions will become available within the 
next few years both from interferometry and from timing. 
The Jodrell Bank data archive will be able to provide values 
or limits on the proper motions of many hundreds of pulsars 
in the near future. These will further improve the constraints 
on the distribution and origin of pulsar velocities. 
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7 CONCLUSIONS 

We have presented an updated catalogue of the kinematics 
of 233 pulsars, effectively doubling the size of the previous 
sample. Our results may be summarised as follows: 

(i) In contrast to the results of Frail, Goss & Whiteoak 
(1994), most pulsars within supernova remnants are found to 
have lower velocities than other pulsars. This is most likely 
due to selection effects where fast moving pulsars leave the 
supernova shell within a relatively short period of time. 

(ii) The proper motions for PSRs B1830-08 and 
B2334+61 are consistent with their proposed associations 
with the supernova remnants W41 and G114.3+0.3. 

(iii) The fastest moving pulsar with a well-defined dis- 
tance is PSR B1133+16 which has a 2-D speed of 640 kms -1 . 
However, according to the CL02 (and TC03) distance model 
PSRs B2011+38 and B2224+65 both have 2-D speeds 
greater than 1500 kms -1 . 

(iv) The CL02 distance model generally predicts smaller 
distances, and hence 2-D speeds, than the TC03 model. The 
mean 1-D and 2-D speeds for pulsars with characteristic 
ages less than 3 Myr are 192(20) and 307(47) kms -1 . The 
observed 1-D and 2-D speeds clearly demonstrate that the 
3-D velocity vector is isotropic. 

(v) Based on a deconvolution analysis of the new sam- 
ples of 1-D and 2-D speeds of young pulsars, we find the 
mean 3-D birth speed to be 400(40) kms - . The 3-D speeds 
are well fit by a Maxwellian distribution with 1-D rms 
a — 265 km s -1 . We find no evidence for a bimodal velocity 
distribution. 
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